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رابطه‌سازی نوین برای جرم ساختگی در روش رهابی پویای لزج * 


ایمان زردی" جواد علامتیان"۲ 


چکیده در اين مفاله رابطه‌سازی نوینی برای جرم ساخنگی در روش رهایی پویای لزج پیشنهاد می‌گردد. نحست معادله‌های حطا برای 
تکرارهای رهایی پویا ارائه می‌شود. سپس, بااستفاده از نگر گر شگورین بهبود يافته, کران‌های جدیدی برای کمیت جرم ساختنگی به دست 
می‌آید. رابطه سازی پیشنهادی به الگوریتم جدیدی برای روش رهایی پویای لزج منجر می‌گردد. برای سنجش کارایی عددی فرایند پیشنهادی» 
سازه‌های محتلفی مانند حرپا و قاب‌های دو بعدی و سه‌بعدی, با رفتارهای حطی و غیرنعطی هندسی تحلیل می‌شوند. نتایج این تحلیل‌ها نشان 
می‌دهد که استفاده از الگوریتم ارائه‌شده و رابطةٌ نوین پیشنهادی برای جرم سانگی, نرخ همگرایی روش رهایی پویای لزج را بهبود می‌دهد. 
به گونه‌ای که همکرایی با تعداد کمتری تکرار انجام می‌پذیرد. 


واژه‌های کلیدی روش رهایی پویاء میرایی لزج. جرم ساختگی, رفتارغیرخطی نیروی نامیزان نرخ همگرایی. 
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مقد مه 

هدف از تحلیل یک سازه. محاسبةٌ تغییرمکان‌ها و تنش- 
های پدید آمده در عضوهای آن, تحت اثر بارهای اعمالی 
تاک و انا ال سا ها را ریت 
تغییرمکان‌های گرهی نوشت. درنتیجه. تغییرمکان‌های 
گرهی به‌عنوان مجهول اصلی از تحلیل سازه به‌دست می- 
آیند. برای رسیدن به این هدف. دو نوع تحلیل خطی و 
غیرخطی (غیرخطی مواد. هندسی و تکیه‌گاهی) وجود 
دارد. اگر ماتریس سختی مستقل از تغییرمکان‌های گرهی 
باشد. تحلیل خطی است. در غیر این‌صورت. تحلیل 
غیرخطی می‌باشد. باید دانست. چنانچه ماتریس سختی 
سازه در هنگام تحلیل تغییر کند. رفتار آن سازه غیرخطی 
خواهد بود. منشاً ایجاد رفتارهای غیرخطی را می‌توان به 
سه دستهٌ غیرخطی مواد هندسی و تکیه‌گاهی تقسیم کرد. 
عامل غیرخطی مواد ناشی از ویژگی‌های ذاتی مصالح 
سازه و تغییرشکل‌های مومسان آن است. رابطةٌ غیررخحطی 
تنش- کرنش, جاری‌شدن مواد و تغییرشکل‌های ماندگار 
ناشی از باربرداری سبب ایجاد این رفتار می‌گردند. 
ازسوی‌دیگر» در حالت غیرخطی هندسی تغییرشکل‌های 
سازه بزرگ درنظر گرفته می‌شود. در این نوع تحلیل 
شکل هندسی سازه در هر مرحله تغییر می‌کند. این رفتار 
هنگامی پدید می‌آید که تغییرمکان‌ها. دوران‌ها و یا 
کرنش‌های سازه بزرگ باشند. چنانچه شرایط تکیه‌گاهی 
سازه در هنگام بارگذاری تغییر کننده رفتار غیرخطی 
تکیه‌گاهی پدید می‌آید. در این پژوهش. نمونه‌های عددی 
با رفتار غیر خطی هندسی تحلیل می‌شوند. دستگاه معادله- 
های حاکم بر رفتار سازه‌های ایستا که بااستفاده از رابطه- 
سازی اجزای محدود يا تفاوت‌های محدود به‌دست 
می‌آید. به‌صورت زیر است: 
)۱( ۳-۳ - 50 

ی را رابطه [ 5 |- 5 ماتریس سختی سازه است 
که متقارن و به ابعاد 0*0 می‌باشد. کمیت ۰1 تعداد 
فوتهای» اراد شاف سک هی 0 


()- و |()عط به‌ترئیب» بردارهای تغییرمکان. 


نشرية مهندسی عمران فردوسی 


رابطه‌سازی نوین برای جرم سانعتگی در روش رهایی پویای لزج 


نیروی داخلی و بار خارجی سازه می‌باشند. برای تحلیل 
غیرخطی سازه‌ها. دو رویکرد وجود دارد؛ روش‌های 
صریح و ضمنی. فرایند تکراری رهایی پویا که در دس 
روش‌های صریح قرار دارد. از عملیات برداری برای 
تحلیل سازه و یافتن پاسخ دستگاه معادله‌های (۱) بهره 
انم وخ کین سیب افتانی کاران» ای بانط 
موردنیاز و سادگی محاسبه‌ها می‌گردد. روش رهایی پویا؛ 
براساس وجود يا عدم وجود میرایی به‌ترتیب» در دو 
دستةٌ روش رهایی پویای لزج و جنبشی جا می‌گیرد. در 
روش رهایی پویای جنبشی. سازه به یک محیط دینامیکی 
منجازی (ساختگی » که فاقد.میرایی است».مشقل می‌شود. 
به‌دلیل نبودن نیروهای میرایی در این محیط دینامیکی 
مجازی» فرض پایستار بودن انرژی مکانیکی سازه برقرار 
است؛ زیرا همه نیروهای موجود در این فضای دینامیکی 
مجازی (شامل نیروهای داخلی اعضا و نیروهای اینرسی)؛ 
پایستار می‌باشند. در روش رهایی پویای جنبشی از 
ویژگی پایستاری این سامانه دینامیکی مجازی استفاده 
می‌گردد؛ به‌گونه‌ای که نقاط با بیشینه کارمايةُ جنبشی در 
نزدیکی تعادل ایستا قرار می‌گیرند [1]. در فن رهایی 
پویای لزج با فرض میرایی ساختگی برای سازه. سامانه 
از فضای ایستا به فضای پویای ساختگی منتقل می‌شود. 
به‌دلیل وجود میرایی ساختگی و عدم وجود بار پویا 
پاسخ گذرا پس از مدتی از بین می‌رود و درنهایت پاسخ 
انیت با زار دشتگام باقن ی مان قایرایی با افتال 
سامانه از فضایی ایستا به پویاء رابطة (۱» به رابطة زیر 
تین ی سوو: 
)۲( ۳ ۲۳ - ]9۳+ ترس + ۸۳۳۳ 

در اين رابطه [زس] - ۷۲ و [ع]- ۳ به ترتیب» 
ماتریس‌های جرم و میرایی ساختگی. قطری و به ابعاد 
0 می‌باشند. بالانویس 0 نشان‌دهنده شمارهٌ تکرار است. 
همچنین 1 و 13. به‌ترتیب. بردارهای سرعت و شتاب 
سازه می‌باشند. 

اتیال ۱۳۱۵ که روش ان تا بط ار نها 


دی [3] معرفی گردید. تاکنون از این روش در حل مسائل 


ستال سین ره تک شماره بک: ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


گوناگونی استفاده شده ات پژوهشگران بسیاری این 
روش را گسترش داده‌اند. این پژوهش‌ها را می‌توان در 
جند دسته جا داد. 

در راستای استفاده از روش رهایی پویا در تحلیل 
سازه‌های مختلف, اتر این روش را برای محاسبهُ تنش‌ها 
و تغییرمکان‌ها در یک مخزن تحت فشار به‌کار برد [4]. 
همچنین» وی توانست یک مسئلة سه‌بعدی ارتجاعی- 
پلاستیک سد قوسی را بااستفاده از روش رهایی پویا حل 
کند [2]. برو و بروتون فن رهایی پویا را برای تحلیل 
خطی و غیر حطی قاب‌ها استفاده کردند [5]. آنها با مقايسة 
راهکار رهایی پویا با روش‌های دیگر برتری این روش 
را در تحلیل خطی و غیرخطی قاب‌ها اثبات نمودند. 
لویس و جونز باکمک روش رهایی پویء به بررسی برهم 


کنش بین پوشش و شبکه در سقف‌های کابلی پیش‌تنیده 
پرداتضتد 161 قرو ستال‌های ار کشت و همکاراتشن 


براساس تئوری خزش. روش رهایی پوبای اصلاح‌شده‌ای 
برای تحلیل ایستا و غیرخطی گنبدهای کابلی, ارائه 
نمودند [7]. 

از فن رهایی پویا در تحلیل صفحات و پوسته‌ها 
بهره برده شده است. رشتن از آن, در تحلیل صفحات 
ارتجاعی که تحت بارگذاری جانبی هستند. بهره برد [8]. 
بر این اساس. تغییرشکل‌ها و لنگرهای ایجاد شده در 
صفحه به‌دست می‌آیند. آلوار و راماچاندرا فن رهایی پویا 
را در مسئلهٌ تغییرشکل‌های بزرگ صفحات ایزوتروییک 
به‌کار بردند [9]. رامش و کریشنامورتی با درنظر گرفتن 
تغییرشکل‌های بزرگ. دوران‌ها و کرنش‌های کوچک به 
تحلیل صفحات و پوسته‌ها با بهره‌گیری از روش رهایی 
پویا پرداختند [110. بارنس توانست بااستفاده از روش 
رهایی پویاء فرایندی عددی برای پیش‌بینی رفتار کابل‌های 
با دهانٌ زیاد و همچنین پوسته‌های شبکه‌ای و غشاهای 
پیش ‌تنیده پيشنهاد دهد [111. 

فن رهایی پوبا در مسائل مربوط به کمانش و 
پس‌کمانش سازه‌ها به‌کار رفته است. رامش و 


کریشنامورتی بااستفاده از روش رهایی پویاء الگوریتم 


تاش هک ان نی ۱۳۹۱۲۰ 


جدیدی برای تحلیل پس کمانشی ارتجاعی و غیرارتجاعی 
سازه‌های خرپایی پيشنهاد کردند [112. برای تخمین 
رظان جر ختتفایق+ارتحاشی کت اقب ترزهای 
فراع اس اشفا تانکبی ونان رنطتهای 
۱ 
روش طول قوس و استفاده از فن رهایی پویای جنبشی. 
تالا تداعس سازهها زا پریشی سیف ۳ 
الگوریتم‌های گوناگونی برای روش رهایی پویا ارائه 
شده است. همچنین, رابطه‌های مختلفی برای تخمین 
عامل‌های تکراری موجود در آن پيشنهاد گردیده است. 
بنس بااستفاده از اید خود که استفاده از اصل ریلی و 
بردار نیروی تابع جرم بود. توانست تخمین مناسبی برای 
غاما شوایس رای ارف دهاز ولن کار و هایر از 
توری گرشگورین برای تخمین عامل‌های متغیر در روش 
رهایی پویا بهره بردند [16]. بااستفاده از تحلیل خطا و 
نگرة گرشگورین, پاپادراکاکیس فرایندی خودکار برای 
تخمین عامل‌های متغیر و تکراری روش رهایی پویا ارائه 
نمود [17]. برپایةٌ گسترش دنبالهٌ تیلون پژند و حکاک 
زانظه‌های وان کید کهشیی کاهی ‏ عکرارها دزی وش 
رهایی پویا گردید [18]. توروی و همکارانش با 
کمینه‌سازی نیروی باقی‌مانده پس از هر تکرار گام زمانی 
اصلاح‌شده‌ای ارائه نمودند [19]. رضایی‌پژند و علامتیان 
الگوریتم تصحیح‌شده‌ای براساس کمینه‌سازی خطای 
تغییرمکان بین دو تکرار پیاپی پيشنهاد کردند [20]. 
رضایی‌پژند و سرافرازی الگوریتم جدیدی برای روش 
رهایی پویا پيشنهاد نمودند به‌گونه‌ای که نیازی به ماتریس 
قبرآیی ی عاما‌های تترضت: تذاشت: زا تانفاو از 
فرایند تکراری استودولا و کمینه‌سازی خطا بین دو تکرار 
پیاپی. رضایی‌پژند و همکارانش روشی برای به‌دست 
آوردن جرم و میرایی ساختگی بهینه معرفی نمودند [22]. 
همچنین» رضایی‌پژند و همکارانش با کمینه‌کردن نیروی 
باقی‌مانده. توانستند گام زمانی جدیدی ارائه کنند 
[23]. علامتیان براساس اجرای تحلیل نموی» رابطهٌ 


جدیدی برای جرم ساختگی در الگوریتم روش رهایی 


نشریه مهندسی عمران فردوسی 


پویای جنبشی ارائه نمود که نرخ همگرایی را افزایش 
می‌داد [1] 

در این مقاله. نخست رابطه‌سازی تکراری برای 
دست‌یابی به معادله‌های خطا در روش رهایی پوبای لزج 
ارائه مین گردنك: سپس. با ترکیب نگرة دایره‌های 
گرشگورین بهبودیافته رابطةٌ نوینی برای جرم ساختگی 
در روش رهایی پویای لزج ارائه می‌گردد [1]. بر این 
اساس» الگوریتم جدیدی برای روش رهایی پویای لزج 
معرفی می‌گردد. درپایان. کارایی رابطه‌سازی و الگوریتم 
پیشنهادی, با تحلیل خطی و غیرخطی سازه‌های مختلف 


بررسی می‌شود. 


رابطه‌های روش رهایی پویای لزج 
روش رهایی پویا یک فرایند تابع اولیه‌گیری عددی است 
که برای یافتن پاسخ‌های دستگاه معادله (۱) به‌کار می‌رود. 
در این روش با داشتن داده‌ها در یک گام کمیت‌های 
مجهول در گام بعدی به‌دست می‌آید. تابع اولیه‌گیری 
عددی بااستفاده از روش تفاوت‌های محدود مرکزی 
انجام می‌پذیرد. بر این اساسء شتاب درجه آزادی [ ام 
عبارت است از تغییر سرعت آن درجه آزادی در بازه 
زمانی بین دو نمو. کمیت تغییر سرعت درجه آزادی ذام 


1 
۰ 1-2 
2 


1 
بین دو نمو زمانی به‌صورت 7 4 -2 8 می‌باشد. بازه 
زمانی تغییر سرعت نیز میانگین دو گام زمانی پیاپی یعنی 


یوخ است. با تقسیم این دو کمیت بر یکدیگر» می- 
توان نوشت: 


عمون. وه 2 
۳( 29یا ,(2 2-4 ) 


۳ 

در اینجا. "۲ و "۰2۳ به‌ترتیب. گام‌های زمانی 2 ام 

و (1+)ام می‌باشند. از سوی دیگره می‌توان بردار نیروی 

نامیزان را که برابر تفاوت نیروی خارجی و داخلی است. 
با بهره گیری از رابطه‌های (۱) و (۲) به‌دست آورد: 


(ع 5۳ - ۲۳ - ۳۳ - ترس + ۷۳1۳ 


نشرية مهندسی عمران فردوسی 


رابطه‌سازی نوین برای جرم سانعتگی در روش رهایی پویای لزج 


در اینج |" < 14 ۳و 1 به‌تر تیب بردارهای نیروی 


نامیزان. بار خارجی و نیروی داخلی وارد به سازه و به 
ایعاد 1و می‌باشند. با جایگزینی رابطة (۳) فر () و 


سم 1 
به‌دست آوردن 2 ب» می‌توان نوشت: 


1 1 +1 
۳+۳ 


0+ 


1 2 - ۳ 2+1 6 


12:0 ز 


بااستفاده از فرایند میانگین گیری سرعت‌ها و رابطة 


(۵ می‌توان نوشت: 


1 دك 2 1 
0 ۱ 
تج 2 + 2 : 
1 1 1 1 
0 ...و2 مب < 1 


شایان توجه است. ماتریس جرم در روش رهایی 
پویا قطری فرض شله است. درنتیجه نیازی به محاسبهةً 
وارون آن نیست. این نکته, سبب کاهش محاسبه‌ها می- 
شود. در رابطه‌سازی روش رهایی پویا از فن تفاوت‌های 
محدود مرکزی استفاده می‌شود. براین اساس می‌توان 


1 
و 
02 12 عّز/, 2 بلج + ۳ 


رابطهٌ (۷) را می‌توان به‌صورت برداری برای کل 
سازه نوشت: 
1 
۸ 2 + - 7 
با بهره‌گیری از رابطه‌های (۵) و (۷ نتيجه زیر 


به‌دست ق ال 


1+1 ۳ 1+1 ۱۳ 
0 2 زگ و 1 
[ِ بس) 1 +( +0 :0+۲0 0 
2۳ 1 ۳ 
(٩)‏ 
+1 +1 
و( ۲+ ۲ ۲ 
لارنص ‏ یز 1 ۳۹ 


با دیفرانسیل گیری از رابطةٌ )٩(‏ می‌توان نوشت: 


+ 1+۱ +2 ۳ 
4 ۳ 7 1 
- )00۳1 +(-+000- 3 +00۳ 
۲ 21۳ 1 210 
1+1 رم 1+1 
4 ۲ 
1 ج( ۲ کب 
۰2( 2 


تال سن ره تک شماره یکت: ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


در اینجاء ()4 به‌معنای دیفرانسیل می‌باشد. باید دانست. 
در رابطه‌سازی روش رهایی پوی. دیفرانسیل گیری در بازه 
کلی همگرایی انجام نمی‌شود؛ بلکه اين فرایند در هر نمو 
زمانی انجام می‌پذیرد. در اين راستاء کمیت‌های روش 
رهایی پویا شامل. جرم میرایی و گام زمانی در هر نمو 
ثابت فرض می‌شوند و تنها پاسخ سازه تغییر می‌کند. این 
شیوه رابطه‌سازی به اين دلیل انجام می‌شود که محیط 
دینامیکی در روش رهایی پویا ساختگی می‌باشد. بنابراین؛ 
برای سادگی هرچه بیشتر محاسبات و صریح شدن آنها؛ 
می‌توان انواع فرض‌ها را در این محیط دینامیکی مجازی 
اعمال کرد. این گونه است که فرض‌هایی مانند قطری‌بودن 
ماتریس‌های جرم و میرایی و نیز ثابت بودن کمیت‌ها در 
هر نمو زمانی معنی می‌یابند. به عبارت دیگر. در هر نمو 
از روش رهایی پویاء با ثابت پنداشتن کمیت‌های جرم. 
میرایی و گام زمانی» شرایطی را فراهم می‌آورند تا 
تغییرمکان سازه. در کوتاه‌ترین زمان ممکن, به پاسخ 
حالت پایدار (استاتیکی) همگرا شود. ازآنجاکه تغییرات 
بار خارجی در تکرارهای رهایی پویا صفر است. بنابراین 
0- ,0 و تغییرات نیروی نامیزان برابر است با: 


۰ ۶۳ ۳ « ۰ 
رد6 - 2( 46 - یذ 
از 1 ال[ 
2,0 ,1 ز 
که در آن. » 8 برابر ژُذ امین درایه ماتریس 
0 
سختی مماسی سازه و می‌باشد. با حطی پنداشتن تغییرات 
نموی تغییرمکان بین دو تکرار پیاپی می‌توان نوشت: 
٩۵۳۱ - 003 , 00۳ 100, ۱-12. 1‏ 


کمیت ,1 نسبت نمو تغییرمکان درجه آزادی 1 ام و 
یا همان نمو خطای تغییرمکان است. باید دانست. خطا 
در گام طام به‌صورت 0۳۰ ۲[ تعریف می‌شود که 
ر 1 ی پاسخ درجة آزادی ز ام است. با 
دیفرانسیل گیری از این خطاء نمو خطای تغییرمکان و نمو 
تغییرمکان با یکدیگر برابر خواهند شد؛ زیرا پاسخ 
دستگاه یگانه است و دیفرانسیل آن صفر می‌شود. برای 
تضمین پایداری روش. باید قدرمطلق نمو خطای 


تال شوه نکن ماه یک :۱۳۹۷ 


تغییرمکان از یک کوچک‌تر باشد تا نمو تغییرمکان به صفر 
میل کند و سیستم به وضعیت ایستا برسد. با تعریف 
نسبت گام زمانی ». به‌صورت: 

۱۳( که 


و بافرض قطری بودن ماتریس جرم ساختگی و 
انجام رابطه‌سازی براساس مقدار ویژه» می‌توان نوشت: 
۹ 


آ 
11 


(0۳ 11:2, (۱( 


10 
که در آن م مقدار وی ماتریس ۳(۳۸6) یا همان 
مجذور فرکانس طبیعی سیستم پویای ساختگی می‌باشد 
[1]. با بهره‌گیری از رابطه‌های (۱۶) و (۱۲) در رابطة 
(۱۰» معادلةٌ درجه دو زیر برای خطای تغییرمکان 

به‌دست اب 
(۱۵) . 0- وبا سوب پر 0 
12 
پراساس تعریف آندروود رابطهٌ بین ماتریس میرایی 
و جرم ساختگی در گام تام می‌تواند به‌صورت زیر بیان 
شود [24]: 
۱1( 112 


رن ۲ تنم 


1 1 


برای ساده‌سازی رابطه‌هاء از تغییرمتغیر /استفاده می - 


سود. ۲ 
سك ۰ 
)۱۷( 12:0 سس گت 7 بت نم سل - 8 
21۳0 


در اين رابطه‌هاء "۳ عامل میرایی در گام « ام 
می‌باشد. برای بررسی شرایط تضمین پایداری» نیاز به 
ريشهٌ معادلهٌ (۱۵) که همان خحطای تغییرمکان است. می- 
باشد. با ساده‌سازی به‌کمک رابطه‌های (۱۷) و (۱۷) و 
فرض 21 ۲۳ ریشه معادلة (۱۵) برابر است با: 


/ 0001 _ 
و 0 +01 3 


1 


(۱۸ 0 ...2 م1 < 1 


ى‌ 
+20 
۱ 2 ۲ 


همان‌طور که مشخص است» ۹ تابعی انتنت از بم. 
۳ و 2. برای رسم نمودار تغیبرات 1 می‌توان یکی 
از دو مقدار *» و 2 را ثابت و دیگری را متغیر درنظر 


نشریه مهندسی عمران فردوسی 


گرفت. اگر دلنای معادلة (۱۵) مثبت باشد» دو ريشة 
حقیقی وجود دارد. برای رسم نمودارها از مقدار بزرگ‌تر 
استفاده شده است. دلیل این موضوع این است که کاهش 
خطا در روش رهایی پویا نموی می‌باشد؛ درنتیجه. نمو 
خطای تغییرمکان به‌تدریج کاهش می‌یابد. ازاین‌ری 
هنگامی که دو عامل برای نمو خطای تغییرمکان وجود 
داشته باشد. کمیت بزرگ‌تر حاکم خواهد بود به‌گونه‌ای 
که موقعیت جدید درج؛ آزادی مورد نظر, به مکان قبلی- 
اش نزدیک‌تر باشد. اگر دلتا مساوی صفر و يا کوچکتر 
از صفر باشد. به ترتیب یک ريشه مضاعف و دو ريشةٌ 
مختلط وجود دارد. 

شکل (۱). تغییرات خطای تغییرمکان (1) را 
دربرابر نسبت گام زمانی ». برای 0 "0 نشان می‌دهد. بر 
این اساس» خطای تغییرمکان ,1 زمانی کمینه می‌شود که 
دلتا برابر با صفر شود. در شکل (۱) با کاهش مقدار ویژه 
هر دو کمیت گام زمانی بهینه و حطای تغییرمکان افزایش 
می‌پابند. به‌عنوان نمونه. با افزایش مقدار ویژه از به ۱۰ 
مقدار نمو خطای تغییرمکان از حدود ۰/۶ به ۰/۳ کاهش 
می‌یابد. از سوی دیگره براساس تعریف مقدار ویژه که با 
رابطه (۱۶) انجام شده است. مقدار ویژه و جرم با یکدیگر 
رابطةٌ وارون دارند؛ یعنی با افزایش مقدار ویژه. کمیت 
جرم کاهش می‌یابد. بنابراین, در هر درج آزادی» هرچه 
مقدار ویژه بزرگ‌تر شود. جرم ساختگی آن کمتر می- 
گردد. درنتیجه. خطای تغییرمکان کوچک‌تر می‌شود. این 
نکته اهمیت درنظر گرفتن جرم ساختگی کمینه را برای 
درجه‌های آزادی نشان می‌دهد. 

اگر مقدار ویژه برای درج؛ آزادی #ام متغیر و عامل 
میرایی ساختگی برابر با یک درنظر گرفته شود شکل (۲) 
به‌دست می‌آید. همانند شکل (۱ با کاهش مقدار ویژه 
گام زمانی بهینه افزایش و خطای تغییرمکان بیشتر شده 
است. 

چنانچه مقدار ویژه برای درجة آزادی ام ثابت و 
برابر با یک و عامل میرایی ساختگی متغیر درنظر گرفته شود. 


شکل ۳ در دسترس قرار مین گیرفن در روش رهایی 


نشريهة مهندسی عمران فردوسی 


رابطه‌سازی نوین برای جرم ساختگی در روش رهایی پویای لزج 


- > #۵ > 0۹۱ 


0 02 04 13 08 1 


پٍ ۳ ۳ ۳ 
1 08 06 04 02 


0 
شکل ۳ تغییرات خطای تغییرمکان برای 2 > ی 0 و 
1< 2 


شکل (۲) برای مقدارهای مثبت عامل میرایی رسم 
شده است. براساس شکل (۳ با افزايش عامل میرایی 
ساختگی. نسبت گام زمانی بهینه افزایش و کمترین 
خحطای تغییرمکان کوچک‌تر می‌شود. 


ستال رشن ی بکا: شماوه یکت ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


اکنون می‌توان شرط پایداری را به‌دست آورد. 
مساوی با صفر شود. با به‌دست آوردن مقدار ویژه( ,1) از 


5 قر [ ]+21 +0 +1 ]2 _ 
( +0(ه) 3 
با جای‌گذاری 2 از رابطهٌ (۲۰) در (۱۸» نسبت 


نموی تغییرمکان‌ها به‌دست اک 


ی 


(۲۱ 0 ...و2 م1 < 1 1 
(1+۵) 


شرط لازم برای پایداری تکرارها این است که قدر 
مطلق خطای تغییرمکان همواره از یک کوچک‌تر باشد. 
بر این اساس. محدودةٌ زیر که در آن پایداری روش 
رهایی پویای لزج تضمین شود به‌دست می‌آید: 
(۲۲( +1> 0 

برای به‌دست آوردن مرزهای ,2 از رابطة (۲۰) 
استفاده می‌شود. بنابراین. محدودة قابل قبول برای 
مقدارهای ویرهٌ سیستم پویای ساختگی که در آن پایداری 


روش تضمین می‌شود. به‌صورت زیر است: 


21 0+ (-21+ 0 -0([ 


(0 +1) 
ث" 
(۲۳) ِ نب ۳ کِ 
[ 2 - 8()0 +201 + رو + 20 
هه 1 
ِ (ه +1) 


نگره‌ی دایره‌های گرشگورین بهبودیافته 
کین از راه‌های تخمین مقدارهای ویژه ماتریس‌هاء استفاده 
از نگرة دایره‌های گرشگورین می‌باشد. در روش رهایی 
پویا مقدارهای ویژه. همان دوره‌های تناوب طبیعی 
نشکا دیتانیکی با کی می‌باشید 11 اهاز این 
نگره و درنظر گرفتن شرایط پایداری, می‌توان ماتریس 
جرم ساختگی را به‌دست آورد. این روش دارای پایه‌های 


تال شوه که ماه :۱۳۹۲ 


۳۱ 


ریاضی قوی انتیت و سبب کاهش زمان محاسبات می - 
شود. این نگره مقدارهای ویژه ماتریس را به‌صورت زیر 
ارائه می‌دهد: 


۱ ۹ 1 ۳( 
(۲( 0 و م2 م1 <- 1 25 ِ > 
اح[ 


۳: 


0 ۹ 


۳ 
َزٌ 


2 


1 


نز 

نگرة دایره‌های گرشگورین یکی از برترین روش - 
ماش مقس فا رای وه سای شا ابا ان کی 
دارای نقاط ضعفی نیز می‌باشد. در این راستا می‌توان به 
ارات زار های قشی یا صفر را فزارهای وه 
ماتریس. اشاره کرد. علامتیان توانست با تقسیم نگرة 
دایره‌های گرشگوزین به سه محدودهٌ زیر نقطةهٌ ضعف 


نگرة گرشگورین را برطرف کند [1. 


ی 
6 مخلردی 17 5 2 < 5 
ِ 
۹ 
اه ٩ > ٩‏ ک ارزو 
(۲2( محدودهٌ 11 ؛ 1 ید 1 اه 
اج 1 
۹ ۳ 
۲۷( محدوده الا : بر - زگ یهد > 5 
از آز 


براساس نگرة دایره‌های گرشگورین بهبود یافته. می - 


توان مقدار ویژه ماتریس را به‌صورت زير تخمین زد: 


۱ ۹ ۳( اک ی 
زر نود ازیو ِ > ,> | زگ - 5 لٍ 
از نز 1 


1 ۱ 1 
از 


1: 


۳ كِ (ِ ِ 0 
 )۲۵(‏ محدودة اا: نک 2 زا نتون5 2 4 
اعز اعز 


0 10 
1 1 
از 
0 ... ,2 م1 < 1 
ف 1۱:2 12 
5 ۵( رب و( 
(۳۸۰) . محدود؛ اا: ز أ > رز> رفس 
۳ 5 
احز 11 11 
,12:۰ < ز 


ا زانجاکه لازم است رابطه‌های (۰)۲۸ (۲۹) و (۳۰) 


شرط مرزی (۲۳) را برقرار کنند. محدودة زیر برای جرم 
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۳۲ 


ساختگی به‌دست می‌آید: 


محدوده 1 : 
۲( +1) 7 > و1 0+0 
[ 0+۵۵ 2- ره +0 ]2 ار نان [ 0+۵۸ + مج ]2 
(۳۱( 12:0 - ز ا< ۱ 
محدودة ]]: 
۲( +) ۲( +1) 


سسسسسسس لت > روز > | 5 9 یت 
[ ۵-0 +200 - ره +6 ]2 ۷ ۱ 0 - 0( +2۰0 + ره +0 [2 


_ 


۳۲( :12 
محدود ا1: 
۲( +1) ۲( +1) 
سس ۰۰ رین 
| (2/1+8(6-8-(۱+0) [2 > ( -06 +2 +( + [2 
بترم اک[ 
(۳۳ 


که وس تست تا 
کمینه نیروی نامیزان می‌شود. بنابراین مرزهای پایینی 
رابطه‌های (۰۳۱ (۳۲) و (۳۳) برای تعیین جرم ساختگی 
به‌کار می‌روند. باتوجه به برابری مرزهای پایینی این سه 
رابطه مقدار جرم ساختگی پیشنهادی از رابطهٌ زیر 
یمین 


(۳( 


8 مرن ‌ 0 


ِ ۱ 20 2 


,2ب ع 1 


ازنگ 


برای سنجش پایداری روش رهایی پویاء باید 
2 1 رابطه‌های (۰۳۱ (۳۲) و (۳۳), از مرزهای 


به‌دست 9 


وم12ع 0 > ,۷+ (0()0-0 +1) 0-2 


(۳۵) محدوده ]: 


0-2۲) 0()0-0( +1>0 12 


(۳2) محدوده 11 


1-2۹ 0ک 1 (۲()۵-۵ +10-202+۲(0 


(۳۷) محدودة ]11 : 


در رابطه‌ها؛ کمیت ۷ به‌صورت زیر تعریف می‌شود: 
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رابطه‌سازی نوین برای جرم ساحتگی در روش رهایی پویای لزج 


۳۸ نگ » 


با حل سه معادلهٌ (۰)۳۵ (۳۱) و (۳۷). محدودهة 
نسبت گام زمانی8» به‌دست می‌آید. این نابرابری‌ها؛ 
محدودة پذیرفتنی برای نسبت گام زمانی را که در آن 


پایداری روش تضمین ِِ ارائه می‌دهند: 


/( 4 208-۷ +2 
(۳۹( 
+20 + 28-1۷1 27 ب 
و12 ز 
پرای محدوده 1 
۷/17( +201 - 2077 + 2۷7 +1 


> ۵ > -1+ 27 


1-1۱2: 


(«ع) 1- ی (۵ +20 + 282 + 
برای محدودة 11 
گت 8+ 5+ 4 کِ 7+ 80 + 8 +80 +8 
2 
4+ 8+ 5+ 7 ۶+ 80 + 8 +80 +8 


21۱2۰ [ 
برای محدودهُ 111 
برای افزایش نرخ همگرایی» از کمترین مقدار 
به‌دست‌آمده از رابطه‌های (۸۳۹ (۶۰) و (۶۱ استفاده 
می‌شود: 


۹ (۲) 


2+ 28-1 -20+0( ۱-۷ 
1 


-1+ 2 +20۴ -20+0(۱۰1- 


4+ 8+ 5 + رل (0 +401- 2+ ۷+ 80۷+ (8+ 8+80 
5 نویه /۵ 


1 


محدوده ۱ 


0 < 2 


ده 


| 120 

سبت گام زمانی در روش رهایی پویا از اهمیت 
بالایی برخوردار است؛ زیر این نسبت برای همه درجه- 
های آزادی سازه یکسان است و اثر مستقیمی بر نرخ 
هگا یی رون دزمان یا داز تراسا هر 
ساختگی ( ع) 


نسبت گام زمانی تابعی از عامل میرایی 


تال ش :و یک: شمارة یک» ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


(۸۲), می‌توان جرم ساختگی پیشنهادی را برای هم 
درجه‌های آزادی سازه از رابطةٌ (۳۶) به‌دست آورد. در 
حالت خاص, اگر میرایی ساختگی صفر درنظر گرفته 
شود آنگاه 0<لمی گردد. درنتیجه. رابطةٌ (۳۶) به‌صورت 
زیر تبدیل می‌شود: 


1-12 :۳( 


و دروب ر 


وابر 
دز درم +20 3 


رابطهٌ بالا با عبارت ارائه‌شده در مرجع [1] برای 


1 


جرم ساختگی روش رهایی پویای جنبشی مطابقت دارد. 
بنابراین» رابطه‌سازی‌های انجام‌شده صحیح می‌باشند. 
شایان توجه است. نابرابری (۳۶) رابطه نوینی را برای 
محاسبةٌ جرم ساختگی روش رهایی پویای لزج ارائه می - 
دهد. در اینجاء برای میرایی ساختگی از رابطهٌ زیر استفاده 
می‌شود [20]: 

3 ۳ ۳( ۳ر(0)) 


الگوریتم رهایی پویای لزج 
در بخش قبلی, با بهره‌گیری از نگرة دایره‌های گرشگورین 
بهبودیافته و بررسی پایداری عددی, رابطةٌ نوینی برای 
جرم ساختگی روش رهایی پریای لزج به‌دست آمد 
(رابطة ۳۶). اکنون؛ الگوریتم روش رهایی پویا براساس 
رابطه‌سازی پیشنهادی ارائه می‌گردد. باید دانست. سنجش 
همگرایی در تکراهای رهایی پویا با دو معیار نیروی 
نامیزان و کارمايةُ جنبشی انجام می‌پذیرد. به سخن دیگر 
تکرارها تا آنجایی انجام می‌شوند که یکی از خطاهای 
نیروی نامیزان و یا کارمایه جنبشی, از معیار همگرایی 
کوچک‌تر شود. الگوریتم پیشنهادی برای روش رهایی 

پویای لزج به‌صورت زير است: 

. فرض مقدارهای اولیه برای سرعت (بردار صفر). 
تغییرمکان (بردار صفر يا در صورت موجود بودن؛ 
تغییرمکان همگراشده در نموی قبل). گام زمانی 
ساختگی و نسبت آن (0-2-1) و معیار همگرایی 
برای نیروی نامیزان و کارمایة جنبشی (1.08-6- به 


ال و و نکن ماه :۱۳۹۲۷ 


۳۳ 


و 52108-12). 

کر 0 مه درم انن‌ضووتا یراب ات با کترین 
نسبت گام زمانی (») به‌دست‌آمده از رابطة (۲). در 
غیر این‌صورت 21 می‌باشد. 

۳ ساخت ماتریس سختی مماسی و بردار نیروی داخلی. 

کال یس ها 

۵. محاسبةٌ نیروی نامیزان بااستفاده از رابطةٌ (4). 

7 اگر به > ۱ تحلیل به گام (۱4) می‌رود. در غیر 
اون خررت آدامم دهم ره 

۷ ساخت ماتریس جرم ساختگی از رابطهٌ پیشنهادی 
(۳۶). 

۸ ایجاد ماتریس میرایی ساختگی بااستفاده از رابطةٌ (44). 

٩‏ به‌هنگام‌سازی بردار سرعت ساختگی بااستفاده از رابطة 
10 

۰ من - # 

۱. اگر به > ۱ تحلیل به گام (6) می‌رود 
در غیر این‌صورت ادامه داده می‌شود. 

۱۲ بهنگام‌سازی بردار تغییرمکان بااستفاده از رابطةٌ (۸). 

۳ افزایش شمارندة تکرار 9+1« 

۶ نتایج همین نمو چاپ می‌شود. 

۵. اگر تعداد تکرارها کافی نبود. تحلیل به گام (۱) باز 
می‌گردد. در غیر این‌صورت تحلیل پایان می‌یابد. 

همان طور که مشخص است. در الگوریتم 

پیشنهادی» نسبت گام زمانی بهینه در هرگام به‌صورت 

خودکار انتخاب می گر دد؛ به گونه‌ای که در آغاز هر نمی 

نسبت گام زمانی برابر با یک فرض می‌شود. با این مقدار 

و فرض گام زمانی برابر با یک» نسبت گام زمانی بهینه از 

رابطةٌ (۶۲) به‌دست می‌آید. بنابراین» نسبت گام زمانی در 


الگوریتم پیشنهادی به‌صورت نسبی تغییر می‌کند. 


نمونه‌های عددی 
برای بررسی کارایی راهکار پیشنهادی. نویسندگان برنامة 
رایانه‌ای در نرم‌افزار متلب (1۸1۸3) نوشته‌اند. این 
برنامه براساس الگوریتم پیشنهادی کار می‌کند و از آن 
برای تحلیل خطی و غیرخطی هندسی سازه‌های مختلف 
مانند خرپا و قاب‌های دو و سه‌بعدی در محدودة رفتار 


نشریه مهندسی عمران فردوسی 


۳ 


ارتجاعی استفاده می‌گردد. تحلیل غیر خطی بااستفاده از 
رابطه‌های اجزای محدود و براساس مختصات هم- 
چرخشی انجام می‌گردد [25]. هر یک از نمونه‌ها به سه 
روش تحلیل شده‌اند. تفاوت این روش‌ها در رابطهٌ 
استفاده‌شده برای جرم ساختگی می‌باشد. با اين فرایند. 
کارایی رابطة پیشنهادی برای جرم ساختگی مشخص می - 
شود. در هر سه روش از رابطهٌ (46) برای محاسبه عامل 
میرایی ساختگی در تکرارهای رهایی پویای لزج بهره 
گرفته شده است. 
در روش نخست که همان راهکار پیشنهادی است. 
از رابطهٌ پیشنهادی (۳۶). برای جرم ساختگی استفاده شده 
است. روش پیشنهادی به اختصار. کاناع ( ۷۲001160 
ممتافحهامک1 متحصفصر۲ متتقعطونع0) نامیده می‌شود. در 
دومین روش که ۳:0۲ (۱۵۵60 224 وه 
صمتای‌هام؟ متصمصون. ععتموصع) نام دار د. جرم 
ساختگی براساس رابطة رضایی‌پژند و علامتیان محاسبه 
شین رنه [20]: 
زر وف :5 0 ۸۲ 2 :9 
(4۵) ۳ : 
2,۰,0 ,1 < [ 
در سومین شیوه از رابطةٌ جرم ساختگی 
آندروود استفاده شده است [24]. راهکار آندروود 
به‌عنوان متداول‌ترین روش رهایی پویا 01۴ 
(طمتا هام۴ متصصفصون پوتفمتن0) شناخته می‌شود. 


رابطةً جرم ساختگی در روش به‌صورت زیر است: 
۱ ۰ ۹ ۳( ۳ ‌ 
۳ / *«1.1- 80۰ 
[ ِ 4 َزٌ 


,1,2 < ز 


(1) 


از مقايسة تکرارهای همگرایی این روش با دو شیوة 0191 


و 0011 بررسی ین گرگ 


نمونهة ۱: خریای فنری. این سازه. سازه‌ای تک درجه 


آزادی و ترکیبی از فنر و یک عضو خرپایی است که در 


نشرية مهندسی عمران فردوسی 


رابطه‌سازی نوین برای جرم سانعتگی در روش رهایی پویای لزج 


شکل (4) نشان داده شده است [20]. سختی فنر و صلبیت 
محوری عضو فولادی به‌ترتیب برابرند با؛ 10:51 نیوتن بر 
سانتی متر و 44483985.77 نیوتن. ازآنجاکه رفتار این 
سازه غیرخطی است. از رابطه‌های زیر برای به‌دست 
آوردن نیروی داخلی و سختی مماسی سازه مزبور استفاده 


می‌شود [24]: 
۱ و ظ ور 2 
[ه 0-3 009-] (-) (0 ۸۳609 0.5 < (۳)] 
0 10 
(6۷) 
م هو( گه) +طی+ 
10 
۱ ۳ 
(۸: () [0 طه2-ج 0 -](6 )۸ 25 ٩‏ 
1 ۳۳ 
0 "وزه اد + وله 
3 : 


بار واردشده به این سازه در12 گام به آن اعمال می- 
شود؛ به گونه‌ای که افزایش بار در هر نمو برابر با4.4484 
نیوتن می‌باشد. شکل (۵) نمودار بار- تغییرمکان این سازه 
را نشان می‌دهد. تعداد تکرارهای لازم برای همگرایی در 
حدول )۱ درج شده‌اند. 

جدول (۱) نشان می‌دهد رابطه‌سازی پیشنهادی 
سبب افزایش نرخ همگرایی روش رهایی پویا شده است؛ 
به گونه‌ای که تعداد تکرارهای همگرایی روش پیشنهادی 
۴و کمتر از دیگر الگوریتم‌های متداول (011 و 
می‌باشد. درنتیجه رابطه‌سازی پیشنهادی سبب 
افزایش نرخ همگرایی روش رهایی پویا و کاهش زمان 
محاسبه‌ها شده است. 
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تال سن ره تک شماره یکت: ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


تحلیل روش 


(94 
(۴ 


۹9( 


۷ 


۷۸ 


1 ]6۳1 


۸ 


۳۵ 


جدول ۱ تعداد تکرارهای همگرایی خریای فنری 


۱۱ 


شکل ۵ مسیر ایستایی خرپای فنری 


۷6۵ ۱۷۹ 


نمونة ۲: قاب توگل . مطابق شکل (0) اين قاب از دو 
عضو(چوبی) با مدول الاستیسیته. سطح مقطع و ممان 
اینرسی, به‌ترتیب برابر با "07110 نیوتن بر مترمربع» 
۳ مترمربعم و "003750 متر به توان چهار 


ال شوه نکن ماه :۱۳۹۲۷۰ 


شکل 1 قاب توگل 


0029316 


۰ 


۸ 


بهبود همگرا ۱ ۳ )/( 
مجموع 
تکرارها 
۵ ۱ ۱" ۲" تکرار خصهه- رصم خی خرطه 
#صه #صه 
۱ 21 1۳ 3 9 
۲ ۱۳ | ۱۶ | ۱۶ ۱۱۵ ۸۷2۷ ۹۳۸۹ 
ع ع 3 3 ۷ 


باه شه است [22] بار کلن وارد به ان سازهترایر نا 
و یرت اس کون ده ترعله: ترازو اعسال 
می‌گردد. برای مدل‌سازی اجزای محدود. هر عضو این 
قاب به پنج جزء تقسیم شده است. شکل (۷): مسیر 
ایستایی این سازه را برای تغییرمکان عمودی گره میانی و 
در دو حالت تحلیل خحطی و غیرخطی نشان می‌دهد. در 
جدول (۲). تعداد تکرارهای لازم برای همگرایی سه 
روش آورده شده است. 

داده‌های جدول (۲) نشان می‌دهد روش پیشنهادی 
۵ با تعداد تکرارهای کمتری نسبت به شیوه‌های 
۶ و ۵۳8 به پاسخ همگرا شده است. 
این بهبود. سبب کاهش زمان محاسبه‌ها و هزین تحلیل 


می‌شود. 
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شکل ۷ مسیر ایستایی قاب توگل 


نشریه مهندسی عمران فردوسی 


۳۹ رابطه‌سازی نوین برای جرم ساحتگی در روش رهایی پویای لزج 


جدول ۲ تعداد تکرارهای همگرایی قاب توگل 
شماره نمو بارگذاری مجموع بهبود همگرایی(/) 
تحلیل روش 
۱ ۲ ۳ ۶ ۵ 4 ۷ ۸ ۹ ۱ تکراره اه - ۲اه یه 0۶اه 
۴اه 0۴ 
094 ۱۹۹۲ ۱۸۳۷ ۱۳۸ ۱۷۹۹ ۱۷۹۹ ۱۷۳۹۶ ۱۷۹۳ ۱۷۹۲ ۱۷۹۲ ۱۷۹۱ ۱۱۸-۱۰ 
خحطی ۹(4 اک ۱۳:۹ ۱۳۳۶ ۱۷۹ ۱۳۳۹ ۱۷۳۵ ۱۳۳۳ ۱۷۳۳۳ رقف ۱۳۳۲ ۱:۳۱ ۱۳ ۰ 
مه ۷:۷ ۱۷۱ ۷ ۹6 1۹6 ۹۶ 1 ۳ ۳ ۳ ۱۷۳۲ 
094 ۳۰۹ ۳ ۳۹۵ ۳:۲ ۱:۳۶ ۳۹۹۸ ۳۸۳ ۳۷/۶ 934 2-۱۰۵۸ ۰۳:۳۳ 
غیر 
۹(4 ۱۹:۹ ۳ ۳۰۳ ۱۱ ۳۳۳۹ ۳5۹۲ ۳۹۹۳ ۳۷۸ 2۱۹ 32 ۳۹۰ 29/۳ 
(9 وف ۷۳۳ .۸ ۸5۸ ۹۹ ۱۰۳۹۵ ۱۹۷ ۱1-۰۳ ۱۹۳۳ ۱۹۳۹۲ 2۱۸۹۳۲ 
جدول ۳ تعداد تکرارهای همگرایی قاب ساختمانی 
۱ شماره نمو بارگذاری جموع بهبود همگرایی(/) 
تحلیا روش ِ 
۱ ۲ ۳ ۳ ۳ 4 ۷ 4 ۷۰ تکرارها عصس خطه رطس 0و 
9 0 
۲ص ۰۳۹ و 1-۰-۰۳۹۷ 2-۰۹۹ 2-۸۰۹ ۳ 2-۰۹۳ 2۰۷۹ 2-۰۷ ٩-۰ ٩12-۰۰‏ 
خطی کر 2-۹ 19۳۵ 3 22۹۳۸ ۳۸ 2929۹۳۷ 9۹ 2۲  (/‏ 22۱ ۳ 2:1 ۸۲ 
صیه ۱۳۹۹ ۱۹۳/۲ ۱۹۳۹۹۹ ۱۹۳۹۹ ۱۹۳۲ ۱۹۳۹۵ ۱۹۳۹۳۹ ۱۹۳۹۳۲ ۱۹۳۹۳۲ ۱۹۳۶ 1۱۲ 
۱ 93 3۳ ۱۳۵۹۸۹ ۱۳9۹۳۹۰ ۳۵۳۹ ۱۳:۸۶ ۱ 2۱۳۸ ۱۳:۷ ۳۳۸ ۱۳:۸ ۱۷( 
ِ -- ۳۹۸ ۱۳۱ ۳۳۸۰ ۳۳۵۲۵ ۳۳۳۹۶ ۳۳۵۲ ۳۳۶۰ ۳۳۹۰ ۳۲۸۵۹ ۳۳/۸۳ ۳۳۹۹ ۶:۰ 4 
خطی ون ۱۳۳۷۶ ۱۱۷۵ ۱9۱۸ ۱۱۱-۰۶ ۱۱۳۸ ۱۱۳۹۶ ۱۱۱۹۹ ۱۱۹ ۱۱۳۰۳ ۱۳۰۹۰ ۱۱۱۹۹۶ 
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تغیرشکل جانبی[00] 


شکل ۸ قاب ساختمانی شک ۶4 سس آیستای قفاب سناختماتن 


نمونة ۲ قفاب ساختمانی. در این نمونه. یک قاب هر طبقه از ساختمان وارد می‌شود. نمودار نیرو- 


ساختمانی که در شکل (۸) نشان داده شده است. تحلیل 
می‌گردد [20]. این قاب ساختمانی دارای 1 طبقه است. 
ستون‌های سه طبقَهٌ نخست و دوم این ساختمان از مقاطع 
۵0 و ۰۷718<*35 و کلیةٌ تیرهای این ساختمان از 
مقطع ۷۷16<31 ساخته شده‌اند. به‌عنوان تقریبی از بار 
مرده و زنده. نیروی گسترد یکنواخت ۵۰ کیلوگرم بر 
سانتی‌مت به کف هر طبقه از ساختمان اعمال می‌گردد. 
همانطوو که فر شعل(۸) تشن داده:شده ات تیروی 


جانبی معادل با نیروی برش پایه زلزله. ب‌صورت افقی به 


نشرية مهندسی عمران فردوسی 


تغییرمکان در شکل »)٩(‏ برای تغییرمکان طبقهٌ آخر این 
شده شتا حدول (۳ تعداد تکرارهای تحلیل خطی و 


غیرخحطی این سازه را نشان می‌دهد. 


نمونة ۶ : خرپای دویعدی . برای سنجش بهتر راهکار 
پیشنهادی. یک خرپای دوبعدی که ۲۹ عضو و ۲۰ درجهٌ 
آزادی دارد و در شکل (۱۰) نشان داده شده است. تحلیل 
می‌گردد [26]. برای این سازه ۸-630 نیوتن» 


ستال سین ره نک شمارة یک: ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


4 متر و یروی ۴ برابر 3*10 نیوتن 
ناک ایکا تشن دم مر خله مه سازه اعمال له ابیع, 

جدول (۶) تعداد تکرارهای لازم برای همگرایی 
هر یک از روش‌ها را در دوحالت تحلیل خطی و 
غیرخطی مشخص می‌کند. نتایج جدول (۶) نشان می‌دهد 
همگرایی روش پیشنهادی 08101 بیشتر از دیگر شیوه‌ها 
(08 و 008 می‌باشد به‌گونه‌ای که با تعداد تکرار 


کمتر پاسخ به‌دست می‌آید. 


تمونه ۵ برج خرپایی. شکل (۱۲) ی خرپایی ۰طبقه 
را نشان می‌دهد. بااستفاده از الگوریتم‌های مختلف روش 
رهایی پویای لزج» این سازه در حالت‌های حطی و 
غیر حطی هندسی تحلیل می‌گردد [1]. این سازه دارای 
۰ درجه آزادی و ۵۶۰ عضو است. سطح مقطع و 


۳۷ 


سانتی‌مترمربع و 6.89510 مگاپاسکال می‌باشد. دو نوع 
بار به این سازه وارد می‌شود؛ یکی بار ثقلی 20 کیلونیوتن 
به هم گره‌های بالاترین طبقه و دیگری بار جانبی باد که 
با رابطة زیر به گره‌های هر طبقه از سازه اثر می‌کند: 


42 ۱ ص ‌ 
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5340 5340 


1 < 1, 2, .., 301 

در این رابطه ,7 ارتفاع طبقهٌ مورد نظر نسبت به 

تراز تکیه گاه می‌باشد. این سازه در دو حالت کشسان 
مسیر ایستایی تغییرمکان افقی بالای این سازه را نشان 
می‌دهد. همچنین. تعداد تکرارهای همگرایی این سازه با 
الگوریتم‌های مختلف روش رهایی پویا در جدول (۵) 


تسم شده‌اند. 


۹) 


جدول ۶ تعداد تکرارهای همگرایی خرپای دو بعدی 


شماره نمو بارگذاری 


تحلیل وش 
۱ ۲ ۳ ۶ 0 4 
وم | ۳۹۲۵ | ۳۲۰ | ع۳ | ۳۲ [ ۱ ۳ 
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شکل ۱۰ خرپای دو بعدی 


ال شور و نکن مارهب ۱۳۹۲۷۰ 


بهبود همگرایی(/) 
ت‌ ۴ ۳3اه اوه - نام 
۷ ۸ ۹ ۱ تکرارها رصم رصم 
۳۹ ۱ ۳9۹۹ ۳5۹۹ ۳۹۶ 
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شکل ۱۱ مسیر ایستایی خرپای دو بعدی 


نشریه مهندسی عمران فردوسی 
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شکل ۱۲ برج خرپایی 


رابطه‌سازی نوین برای جرم سانعتگی در روش رهایی پویای لزج 


جدول ۵ تعداد تکرارهای همگرایی برج خرپایی 


شماره نمو بارگذاری 
تحلیل روش 
۱ ۲ ۳ ع ۵ 1 

۹ ۸ ۸۵۵۱۸ | ۸۵2۳۹ | ۸۵۱۰ ] ۸۵۳۵۹۵ | ۸۵۳۸۵ 
خطی ۳4 ف ال وا زا ۵ مش ۵ ۸۱۱۷۸ 
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« صوصه | ۳۲۳۸۲۰ ۳۰٩‏ ] ۳۱۱۲۰ | ۳۱۲۳۵ ] ۳۱۳۵۱ ]۳۱۶۷۱۰ 


باتوجه به شمار تکرارهای جدول (۵) می‌توان دریافت در 
هر دو تحلیل خطی و غیرخطی, نرخ همگرایی روش 
پیشنهادی 28 بیشتر از دیگر راهکارها (*[0010 و 
۶ می‌باشد. به سخن دیگر. راهکار پیشنهادی دارای 
کارایی بسیار مناسب در تحلیل سازه‌های بزرگ با تعداد 
زیاد درجة آزادی می‌باشد. 

از سوی دیگر. درصد بهبود همگرایی روش رهایی 
پویای پیشنهادی (1«1ع) و شیوة 011 درمقایسه با 
فرایند 1 برای همه نمونه‌های عددی در شکل (۱۶) 
رسم شده است. شکل (۱۶) برتری نرخ همگرایی شیوة 
پیشنهادی را برای تحلیل انواع سازه‌ها شامل خرپا و قاب- 
های دو و سه‌بعدی نشان می‌دهد. شایان توجه است. در 
نمونهٌ ۱ (خرپای فنری) که سازه یک درجه آزادی دارد؛ 
نسبت گام زمانی تکرارهای رهایی پویا برابر با یک درنظر 


نشرية مهندسی عمران فردوسی 
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جابه جایی تیپ [617] 
شکل ۱۳ مسیر ایستایی برج خرپایی 
بهبود همگرایی(/) 
مجموع 
۷ ۸ ۹ ۱ تکرارها ۴ هام اوه - انامه 
۴و ۴و 
۸۳۷۹ ۸9۵۳۷۵ ۸۹۳۷۱ ۱+ ۱+ 
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گرفته شده است. در نمونه‌های دیگر» نسبت گام زمانی 
هر تکرار از رابطة پیشنهادی (4۲) به‌دست آمده است. 
استفاده از این رابطه برای سبت گام زمانی» سبب نضمین 
همگرایی و همچنین افزايش چشم‌گیر نرخ همگرایی 
نی که 
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اب توگل قاب ساختماتی خریای دو بعدی خریای قتری 


نوع سازه 


بر خریایی 


شکل ۱۶ درصد بهبود همگرایی در نمونه‌های مختلف در تحلیل غیر خطی 


سال سی و یک شمازة یک» ۱۳۹۷ 


ایمان زردی- جواد علامتیان 


نتیجه گیری 
در این مقاله رابطه‌سازی نوینی برای جرم ساختگی در 
روش رهایی پویای لزج پيشنهاد گردید. براین اساس 
لگزرسیس اف بان از ارت رای 
گرشگورین بهبود یافته به‌دست آمد [1]. برای سنجش 
عددی. سازه‌های مختلف مانند خرپا و قاب‌های دو و 
سه‌بعدی با رفتارهای کشسان خطی و غیرخطی هندسی. 


۳۹ 


شیوة پیشنهادی ۵15 که در آن از رابطةٌ جرم ساختگی 
بهیو دیافته استفاده شده است» بیشعر از روش‌های متداول 
۶ و 708 می‌باشد. به سخن دیگره الگوریتم 
پیشنهادی با تعداد تکرار کمتری نسبت به دیگر روش‌ها 
همگرا می‌شود. درنتیجه, زمان محاسبه‌ها و حجم عملیات 
کامپیوتری در فرایند پیشنهادی کمتر از دیگر الگوریتم- 
های رهایی پویای لزج می‌باشد. 


تحلیل شدند. نتایج عددی نشان داد نرخ همگرایی 
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